Structured Error Recovery for Codeword-Stabilized Quantum Codes 
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Codeword stabilized (CWS) codes are, in general, non-additive quantum codes that can correct 
errors by an exhaustive search of different error patterns, similar to the way that we decode classical 
non-linear codes. For an n-qubit quantum code correcting errors on up to t qubits, this brute-force 
approach consecutively tests different errors of weight t or less, and employs a separate n-qubit 
measurement in each test. In this paper, we suggest an error grouping technique that allows to 
simultaneously test large groups of errors in a single measurement. This structured error recovery 
technique exponentially reduces the number of measurements by about 3* times. While it still 
leaves exponentially many measurements for a generic CWS code, the technique is equivalent to 
syndrome-based recovery for the special case of additive CWS codes. 
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I. INTRODUCTION 

Quantum computation makes it possible to achieve 
polynomial complexity for many classical problems that 
are believed to be hard ||. To preserve coherence, 
quantum operations need to be protected by quantum 
error correcting codes (QECCs) ||-[|] . With error proba- 
bilities in elementary gates below a certain threshold, one 
can use multiple layers of encoding (concatenation) to re- 
duce errors at each level and ultimately make arbitrarily- 
long quantum computation possible |6|-|l4| ■ 

The actual value of the threshold error probability 
strongly depends on the assumptions of the error model 
and on the chosen architecture, and presently varies from 
10 _3 % for a chain of qubits with nearest-neighbor cou- 
plings [|l5| and 0.7% for qubits with nearest-neighbor cou- 
plings in two dimensions to 3% with postselection 
jFJj, or even above 10% if additional constraints on errors 
are imposed [ jl3| . 

The quoted estimates have been made using stabilizer 
codes, an important class of codes which originate from 
additive quaternary codes, and have aparticularly simple 
structure based on Abelian groups [n|]. Recently, a 
more general class of codeword stabilizedJCWS) quan- 
tum codes was introduced in Refs. ]l8|-|2l||. This class 
includes stabilizer codes, but is more directly related to 
non-linear classical codes. 

This direct relation to classical codes is, arguably, 
the most important advantage of the CWS framework. 
Specifically, the classical code associated with a given 
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CWS quantum code has to correct certain error patterns 
induced by a graph associated with the code. The graph 
also determines the graph state |22] serving as a starting 
point for an encodingalgorithm exploiting the structure 



of the classical code 20 
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With the help of powerful tech- 
niques from the theory of classical codes, already several 
new families of non-additive codes have been discovered, 
including codes with parameters proven to be superior to 
any stabilizer code ||-|^. 

Both classical additive codes and additive quantum 
codes can be corrected by first finding the syndrome of 
a corrupted vector or quantum state, respectively, and 
then looking up the corresponding error (coset leader) in 
a precomputed table |27j]. This is not the case for non- 
linear codes. In fact, even the notions of a syndrome and 
a coset become invalid for general non-linear codes. Fur- 
thermore, since quantum error correction must preserve 
the original quantum state in all intermediate measure- 
ments, it is more restrictive than many classical algo- 
rithms. Therefore, the design of a useful CWS code must 
be complemented by an efficient quantum error correc- 
tion algorithm. 

The goal of this work is to address this important unre- 
solved problem for binary CWS codes. First, we design 
a procedure to detect an error in a narrower class, the 
union stabilizer (USt) codes, which possess some partial 
group structure [^3[ |3| |2f|. Then, for a general CWS 
code and a set of graph-induced maps of correctable er- 
rors forming a group, we construct an auxiliary USt code 
which is the union of the images of the original CWS 
code shifted by all the elements of the group. Finally, 
we construct Abelian groups associated with correctable 
errors located on certain index sets of qubits. The actual 
error is found by first applying error-detection to locate 
the index set with the relevant auxiliary USt code, then 
using a collection of smaller USt codes to pinpoint the 
error in the group. Since we process large groups of er- 
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rors simultaneously, we make a significant reduction of 
the number of measurements compared with the brute 
force error correction for non-linear (quantum or classi- 
cal) codes. 

More precisely, we consider an arbitrary distance-<i 
CWS code ((n,K,d)) that uses n qubits to encode a 
Hilbert space of dimension K and can correct all t-qubit 
errors, where t = [(d - l)/2j. In Sec. £| we give a brief 
overview of the notations and relevant facts from the the- 
ory of quantum error correction. Then in Sec. [III| , we 
construct a reference recovery algorithm that deals with 
errors individually. This algorithm requires up to B(n, t) 
measurements, where 



B(n,t) 



t 

E 

i=0 



(1) 



is the total number of errors of size up to t (this bound 
is tight for non-degenerate codes). Each of these mea- 
surements requires up to n 2 + KO(n) two-qubit gates. 
In order to eventually reduce the overall complexity, we 
consider the special case of USt codes in Sec. [V. Here 
we design an error-detecting measurement for a USt code 
with a translation set of size K that requires 0(Kn 2 ) 
two-qubit gates to identify a single error. Our error 
grouping technique presented in Sec. [v| utilizes such a 
measurement to check for several errors at once. For 
additive CWS codes the t echnique reduces to stabilizer- 
based reco very [Sec. V_C[. In the case of generic CWS 
codes [Sec. (VD|], we can simultaneously check for all er- 
rors located on size-i qubit clusters; graph-induced maps 
of these errors form groups of size up to 2 2t . Searching 
for errors in blocks of this size requires up to ("J — 1 
measurements to locate the cluster, plus up to 2t addi- 
tional measurements to locate the error inside the group. 
In Sec. [v] we discuss the obtained results and outline 
the directions of further study. Finally, in Appendix ^ 
we consider some details of the structure of corrupted 
spaces for the codes discussed in this work. 

Note that some of the reported results have been pre- 
viously announced in Ref. p9|. 



II. BACKGROUND 
A. Notations 

Throughout the paper, "H 2 = C 2 denotes the complex 
Hilbert space that consists of all possible states a|0)+/3|l) 
of a single qubit, where a,/3 G C and |a| + |/3| = 1. 
Correspondingly, we use the space Uf n = (C 2 )®" = C 2 " 
to represent any n-qubit state. Also, 



1 {I,X,Y,Z} 



0,...,3 



(2) 



denotes the Pauli group of size 2 2n+2 , where X, Y, Z are 
the usual (Hermitian) Pauli matrices and I is the iden- 
tity matrix. The members of this group are called Pauli 



operators; the operators in (|J) with m = form a basis 
of the vector space that consists of all operators acting 
on n-qubit states. The weight wgt(-EJ) of a Pauli operator 
E is the number of terms in the tensor product (|2j) which 
are not a scalar multiple of identity. Up to an overall 
phase, a Pauli operator can be specified in terms of two 
binary strings, v and u, 



U = Z V X U 



Z^Z" 2 



ZV n VUl y1t2 
n A l A 2 



Hermitian operators in £P n have eigenvalues equal to 1 
or — 1. Generally, unitary operators (which can be out- 
side of the Pauli group) which are also Hermitian, i.e., all 
eigenvalues are ±1, will be particularly important in the 
discussion of measurements. We will call these operators 
measurement operators. Indeed, for such an operator M, 
a measurement gives a Boolean outcome and can be con- 
structed with the help of a single ancilla, two Hadamard 
gates, and a controlled M gate ]l6| (see Fig. ||). The al- 
gebra of measurement operators is related to the algebra 
of projection operators discussed in j30|, but the former 
operators, being unitary, are more convenient in circuits. 
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FIG. 1: Measurement of an observable M with all eigenval- 
ues ±1. The first Hadamard gate prepares the ancilla in 
the state (|0) + \ip) /y/2, hence \a) = (|0) + |1)) \ip) j~j2. 
The controlled- M gate returns |6) = C M \a) = (|0) + 
|1) M \tp))/y/2. The second Hadamard gate finishes the in- 
complete measurement, \c) — |0) Pm \4>) + |1) Qm \ip), where 
we used the projector identities (^) . If the outcome of the an- 
cilla measurement is |0), the result is the projection of the ini- 
tial n-qubit state \ip) onto the +1 eigenspace of M (Pm | ■*/>}), 
otherwise it is the projection onto the —1 eigenspace of M 
(Qm IV'})- For an input state |1) \tf>) with ancilla in the state 
|1), the circuit returns |1) P M |V>) + |0) Qm |V>)- 

A measurement of an observable defined by a Pauli op- 
erator M will be also called Pauli measurement |nj . For 
lack of a better term, other measurements will be called 
non-Pauli; typically the corresponding circuits are much 
more complicated than those for Pauli measurements. 

We say that a state \ip) € Hf n is stabilized (anti- 
stabilized) by a measurement operator M if M\ip) = \ip) 
(M\ip) = — The corresponding projectors onto the 
positive and negative eigenspace are denoted by Pm and 
Qm, respectively; they satisfy the identities 



M = P M - Qm = 2P M - 1 = 1 - 2Q 



M ■ 



(3) 



We say that a space Q is stabilized by a set of operators 
M. if each vector in Q is stabilized by each operator in M. . 
We use V(A4) to denote the maximum space stabilized by 
A4, and V ± (A4) to denote the corresponding orthogonal 
complement. For a set Ai of measurement operators, 
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each state in V (M) is anti-stabilized by some operator 
in M. 

When discussing complexity, we will quote the two- 
qubit complexity which just counts the total number of 
two-qubit gates. Thus, we ignore any communication 
overhead, as well as any overhead associated with single- 
qubit gates. For example, the complexity of the mea- 
surement in Fig. |l| is just that of the controlled- M gate 



operating on n + 1 qubits |39|. For all circuits we discuss, 



the total number of gates (single- and two-qubit) is of the 
same order in n as the two-qubit complexity. 



B. General QECCs 

A general n-qubit quantum code Q encoding K quan- 
tum states is a i^-dimensional subspace of the Hilbert 
space Wf n . Let {li)}^ be an orthonormal basis of the 
if-dimensional code Q and let £ C £P n be some set of 
Pauli errors. The overall phase of an error [i m in Eq. ^] 
is irrelevant and will be largely ignored. The code detects 
all errors E € £ if and only if B nffl 



(j\E\i) =C E 6 ij 



(4) 



where Ce only depends on the error E, but is indepen- 
dent of the basis vectors 1 32 . The code has distance d if 
it can detect all Pauli errors of weight (d— 1), but not all 
errors of weight d. Such a code is denoted by ((n, K, d)). 

The necessary and sufficient condition for correcting 
errors in £ is that all non-trivial combinations of errors 
from £ are detectable. This gives ||, H 



{j\E\E 2 \i)=C EuE ,5 ij , 



(5) 



where E±,E 2 £ £ and, again, CelE 2 is the same for all 
basis states i, j. A distance-d code corrects all errors of 
weight s such that 2s < d— 1, that is, s < t = [(d— 1)/2J . 

The code is non-degenerate if linearly independent er- 
rors from £ produce corrupted spaces E(Q) = {E : 
€ 2} whose intersection is trivial (equals to {0}); 
otherwise the code is degenerate [[l7|. A stricter condi- 
tion that the code is pure (with respect to £ ) requires 
that the corrupted spaces Ei(Q) and E 2 {Q) be mutually 
orthogonal for all linearly independent correctable errors 
Ei, E 2 G £ ■ 

For a degenerate code, we call a pair of correctable 
errors E\,E 2 € £ mutually- degenerate if the corrupted 
spaces Ei(Q) and E 2 (Q) coincide. Such errors belong to 
the same degeneracy class. For recovery, one only needs 
to identify the degeneracy class of the error that hap- 
pened. The operators like e\e 2 , connecting mutually- 
degenerate correctable errors Ei and E 2 , have no effect 
on the code and can be ignored. 

As shown in Appendix [a[ for all codes discussed in this 
work, any two correctable errors E\, E 2 yield corrupted 
spaces Ei(Q), E 2 (Q) that are either identical or orthog- 
onal. Then, errors from different degeneracy classes take 



the code to corrupted spaces that are mutually orthog- 
onal. Also, for these codes, a non-degenerate code is al- 
ways pure. In terms of the error correction condition (||) , 
we have Ce ± ,e 2 = for errors E\, E 2 in different de- 
generacy classes and Ce x ,e 2 ^ f° r errors in the same 
degeneracy class. 



C. Stabilizer codes 

Stabilizer codes |l6| are a well known family of quan- 
tum error-correcting codes that are analogous to clas- 
sical linear codes. An [[n, k, d}] stabilizer code maps a 
2 fc -dimensional fc-qubit state space into a 2 fc -dimensional 
subspace of an n-qubit state space. 

The code is defined as the space stabilized by an 
Abelian subgroup of the n-qubit Pauli group, y C & n , 
with n — k Hermitian generators, y — {G\, . . . , G n -k)- 
For such a space to exist, it is necessary that — 1 ^ y. 
The Abelian group y is called the stabilizer of Q. Ex- 
plicitly, 



Q = {ty):S\rl>) = \il>), VSey}. 



(0) 



The code Q is stabilized by y iff it is stabilized by all 
n — k generators G{. In other words, it is an intersection 
of subspaces stabilized by Gi, 



Q = p| V(G t 



(7) 



The normalizer of y in denoted as is the 
group of all Pauli operators U which fix y under con- 
jugation {WSU = S for all S G y). The term nor- 
malizer reflects the fact that these operators commute 
with y Jul) . It is possible to construct 2k logical op- 
erators Xj, Zj, j = I, . . . , k, with the usual commuta- 
tion relations, that together with the generators of y 
generate the normalizer (modulo an overall phase fac- 
tor) |l6|, 53] . The (n — k) generators G t of the stabi- 
lizer, along with the k operators Zj, generate a subgroup 
y = (Gi, . . . , G„_fc, Zi, . . . Zk) of which becomes a 
maximal Abelian subgroup when including the generator 
i%. The group y stabilizes a unique state 



|*> = |0...0>. 



(8) 



The operators Xj acting on \s) generate the basis of the 
code, 



(9) 



Generally, each detectable Pauli error Ej € that 
acts non-trivially on the code anti-commutes with at least 
one generator Gi, and errors from different degeneracy 
classes anti-commute with different subsets of y. We 
can thus identify a degeneracy class by the set of genera- 
tors Gi which anti-commute with it. The corrupted code 
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space Ej(Q) = {Ej\tp) : e Q} is anti-stabilized by 
those generators Gi that anti-commute with Ej. Indeed, 

Gi{Ej\4>)) = (—EjGi)\ij)) = —Ej\ip), 

which means that the measurement Gi of gives — 1. By 
measuring all Gi, we get the syndrome that consists of 
n—k numbers 1 or — 1. There are in total 2™~ fc possible 
syndromes identifying different error degeneracy classes, 
including the trivial error 1 which corresponds to the 
all-one syndrome vector. 

Any two corrupted spaces Ei{Q) and Ej(Q) are mutu- 
ally orthogonal or identical. The whole 2™-dimensional 
n-qubit state space Hf n is thus divided into L = 2 n ~ k 
orthogonal 2 fc -dimensional subspaces Qj =Ej(Q), 

L-l 

Wf* = 0ei, Qi±Qj for i^j. (10) 
j=o 

The representatives of different error classes can be cho- 
sen to commute with each other and with the logical 
operations JQ. These representatives form an Abelian 
group [Q & = (gx, . . . g n -k) whose generators gi can be 
chosen to anti-commute with only one of the generators 
of the stabilizer each, giGj = (— l) 5iJ Gjgi (this follows 
from Proposition 10.4 in Ref. g). Altogether, the gen- 
erators {gx, . . . , g n -k} can be regarded as a set of Pauli 
operators forming the basis of the cosets of the normal- 
izer JV of the code Q in 3P n . 

Example 1. The [[5,1,3]] stabilizer code is defined by 
the generators 

Gx =XZZXI, G 2 = IXZZX, 

(11) 

G 3 =XIXZZ, d = ZXIXZ. v ' 

For this code, the logical operators can be taken as 

X = ZZZZZ, Z = XXXXX. (12) 
A basis of the code space is (up to normalization) 

4 

|0) = II (1 + G,)|00000), \1)=X\0). 

4=1 

By construction, both basis states are stabilized by the 
generators Gi. The corresponding stabilizer group is 
3* = (Gx, ■ ■ ■ ,G^). The group of equivalence classes 
of correctable errors is generated by the representatives 
(note the mixed notation, e.g., ZxZ^ = ZIZII) 

.9i = ZiZz, g 2 = ZZZZI, g 3 = ZZIZZ, g A = Z 2 Z b . 

The gj are chosen to commute with the logical operators 
and also to satisfy Gtgj = (— l) Sij gjGi. Note that the 
operators of weight one forming the correctable error set 
do not by themselves form a group. The generators gi can 
be used to map correctable errors to the corresponding 
group elements with the same syndrome. This gives, e.g., 
Z 2 -> .9254 = Z 2 X, X 2 .91 = ZxZ 3 , Y 2 -> g\g 2 g± = 
Z\Z 2 Z^,X = Z/±Z^. □ 



D. Union Stabilizer codes 

The decomposition ( ^0[ ) can be viewed as a construc- 
tive definition of the Abelian group = (gx, ■ ■ ■ Qn— k) °f 
all translations of the original stabilizer code Q in 3P n . 
(In the following, this stabilizer code is denoted Qq.) Any 
two non-equivalent translations t j ,tj € & belong to dif- 
ferent cosets of the normalizer jV§ of the code Qq in & n 
and, therefore, yield mutually orthogonal shifts 

U(Qo)±t j (Q ), i^j. (14) 

A union stabilizer code |23|, ^ (USt) is a direct sum 

K 

C = 0*i(Go) ( 15 ) 

8=1 

of K shifts of the code Qo by non-equivalent translations 

u g sr, i = i,...K. 

The basis of the code defined by (|l^) is the union of the 
sets of basis vectors of all ti(Qo). As a result, the dimen- 
sion of Q is K 2 k . This code is then denoted ((n,K2 k ,d)), 
where d is the distance of the new code. Generally, this 
distance does not exceed the distance of the original code 
with respect to the same error set, d(Q) < d(Qo). How- 
ever, if the code Q is degenerate and the original code 
Qo is one-dimensional, this need not be true. 

E. Graph states 

The unique state \s) defined by Eq. (|J) itself forms 
a single-state stabilizer code \[n, 0,c?']]. Its stabilizer 
3 s = (Gx, ■ ■ ■ , G n -k, Zx, ■ ■ ■ Zk) has exactly n mutually 
commuting generators. Note that for stabilizer states we 
follow the convention that the distance d' is given by the 
minimum weight of the non-trivial elements of the stabi- 
lizer 3 '. Generally, such a state stabilized by an Abelian 
subgroup in of order 2"_( which does not include —1) 
is called a stabilizer state |lq] . 

A graph state [§2| is a stabilizer state of a group whose 
n mutually commuting generators Si are defined by a 
(simple) n- vertex graph Q with an adjacency matrix R £ 
{0, i}™ xn . Specifically, the generators are 

Si = X t Z^ 2 . . . Zf^ = X t Z r > , (16) 

where r^, i = 1, . . . n denotes the ith row of the adjacency 
matrix R. A graph state can be efficiently generated 
HD by first initializing every qubit in the state |+) = 
(|0) + |1))/V2 (e.g., applying the Hadamard gate on |0)), 
and then using a controlled- Z gate Cfj on every pair 
(i,j) of qubits corresponding to an edge of the graph 
Q = {V,E), 

\ s )= II CZ j H® n \0} n = Ug\0) n , (17) 
where |0)„ is a state with all n qubits in state |0). 
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Any stabilizer state is locally Clifford-equivalent (LC- 
equivalent) to a graph state [j36|-|38f . That is, any sta- 
bilizer state can be transformed into a graph state by 
individual discrete rotations of the qubits. This defines 
the canonical form of a stabilizer state. 





FIG. 2: Ring graphs with (a) 3 vertices and (b) 5 vertices. 

Example 2. Consider the ring graph for n = 3 
[Fig. ^(a)] which defines the stabilizer generators Si — 
XZZ, S 2 = ZXZ, S 3 = ZZX. The corresponding sta- 
bilizer state \s) is an equal superposition of all 2 3 states 
[result of the Hadamard gates in Eq. (|l7|)1, taken with 
positive or negative signs depending on the number of 
pairs of ones at positions connected by the edges of the 
graph [result of the gates Cfj = (— l) y ]. In the following 
expressions we omitted normalization for clarity: 



|000) + |001) 
1 100) - 1 101) 
|010) - |011) 
|110> - Jill) 



|010) - |011) 
|110)-|111) 
|000) + |001) 
1 100) - |101). 



(18) 
(19) 



CWS codes and their standard form 



Codeword stabilized (CWS) codes |1^, |2(j represent a 
general class of non-additive quantum codes that also in- 
cludes stabilizer codes. They can be viewed as USt codes 
originating from a stabilizer state. Any CWS code is lo- 
cally Clifford-equivalent to a CWS code which originates 
from a graph state |2C|| . 

Specifically, a CWS code {{n,K,d)) in standard form 
p0[ is defined by a graph Q with n vertices and a classi- 
cal code C containing K binary codewords Cj. The orig- 
inating stabilizer state is the graph state \s) defined by 
Q, and the codeword operators [translations in Eq. (|l5|)] 
have the form Wi = Z Ci , i = 1,...K. For a CWS code 
in standard form we use the notation Q = (Q,C). 

An important observation made in Ref. pOfl is that a 
single qubit error X, Z , or Y acting on a code state 



Wi\ 



(20) 



is equivalent (up to a phase) to an error composed only 
of Z operators. This establishes the following mapping 
between multi-qubit errors and classical binary errors. 
First, let Ei = Z^ ' 1 ^ Xi be an error acting on the ith 
qubit of Wjls). Then we see that 

EiWj\s) = EiWjSils) = ±(^Sj)W^», 



where the term E;Si 



z {o,i} ZI 



consists only of opera- 



tors Z. The general mapping of an error E 



from the error set £ C 
{0, 1}™ is defined as 



n Z v X u 

3?ti to a classical error vector in 



C\g(E) =v®0 Ui r ; , 



(21) 



i=i 



where ui is the I th component of the vector u. We will 
refer to both the binary vector C\g(E) and the operator 
ZC\ g (E) ag graph i ma g e f the Pauli operator E. 

Theorem 3 from Ref. |2(| establishes the correspon- 
dence between the error-correcting properties of a quan- 
tum code Q and those of the corresponding classical code 
C. It states that a CWS code (in standard form) given by 
a graph Q and word operators {Wi = Z c } ce c detects er- 
rors in the set £ if and only if the classical code C detects 
errors from the set C\g(£), and for each E E £, 



either Clg(E) ^ 0, 
or, for each i, Z C 'E = EZ C 



(22) 
(23) 



The code Q is non-degenerate (and also pure, see Ap- 
pendix |a|) iff condition ( p2|) is satisfied for all errors in £ 
p0| , pi] ]. For a degenerate code, condition (|23l) n eeds to 
be ensured for errors that do not satisfy Eq. (|22|). 

The beauty of the CWS construction is that, for a 
given code in standard form, we no longer need to worry 
about possible degeneracies. The classical error patterns 
induced by the function Clg(-) also separate the errors 
into corresponding degeneracy classes [g(| ^l| . 

The CWS framework is general enough to also include 
all stabilizer codes Specifically, a stabilizer code 

[[n, fc, d]} with the stabilizer y = (Gi, . . . G n -k) and log- 



ical operators Zi, 
the stabilizer 



Xi, corresponds to a CWS code with 



■y = (Gi, 



■ G, 



Z, 



,Z k ) 



(24) 



and the codeword operator set W — (Xi, . . . , Xk) form- 
ing a group of size K — 2 k . Generally, an LC transfor- 
mation is required to obtain standard form of this code. 

Conversely, an additive CWS code Q = ((n,K,d)) 
where the codeword operators form an Abelian group 
(in which case K = 2 k wit h integer k) is a stabilizer code 
[[n, k, d]} ^(J. In Sec. VB we show that the n — k genera- 



tors Gi of the stabilizer can be obtained from the graph- 
state generators Si [Eq. (|l6|)1 by a symplectic Gram- 
Schmidt orthogonahzation procedure p3[ which has no 
effect on the codeword operators. 

Example 3. Consider a non-degenerate CWS code 
((5, 6,2)). The corresponding ring graph is illustrated in 
Fig. |(b) (Ref. @). The n = 5 generators of the sta- 
bilizer .y are S\ = XZIIZ and its four cyclic permuta- 
tions. The corresponding stabilizer state has a structure 
similar to Eq. (|l8|), but with more terms. Word operators 
have the form Z Ci , with the classical codewords 



c = 00000, ci = 01101, c 2 = 10110, 
c 3 = 01011, c 4 = 10101, c 5 = 11010. 



(25) 
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□ 

Example 4. To express the [[5,1,3]] stabilizer code 
[see Example 1] as a ((5,2,3)) CWS code in standard 
form, we explicitly construct alternative generators Si 
of the stabilizer & = (G u G 2 , G 3 , G A ,Z) [Eq. @] to 
contain only one X operator each. We obtain S3 = 
G1G2Z = IZXZI and its four cyclic permutations. 
This does not require any qubit rotations due to a 
slightly unconventional choice of the logical operators in 
Eq. (jl2|). The corresponding graph is the ring jpO| , sec 
Fig. 0(b). The codeword operators are W = {7, A}, 
which by Eq. ( |l^ ) correspond to classical binary code- 
words {00000, 11111}. Note that the error map induced 
by the graph is different from the mapping to group 
elements in Example 1; in particular, Z cle ^ z ' 2 ' ) = Z 2 , 

z c\ g (x 2 ) = z t z 3 , z a e<- Y ^ = z 1 z 2 z 3 . a 



III. GENERIC RECOVERY FOR CWS CODES 

In this section we construct a generic recovery algo- 
rithm which can be adapted to any non-additive code. 
To our best knowledge, such an algorithm has not been 
explicitly constructed before. 

The basic idea is to construct a non-Pauli measurement 
operator Mq — 2Pq — 1, where Pq = J2iLi \ w i)( w i\ is 
the projector onto the code Q spanned by the orthonor- 
mal basis {\wi)}fL 1 . The measurement operator is fur- 
ther decomposed using the identity 



K 



K 



-M Q = 1 -2^K>{ttj| - ]J(t - 2\ Wi ){wi\). (26) 

i i 

We use the graph state encoding unitary Ug from 
Eq. (|l7]) and the following decomposition [|o| of the 
standard-form basis of the CWS code in terms of the 
graph state \s) and the classical states \ci): 



Z c >\s) 



Z c *Ug\0) n = UgX c *\0) n = Ug\ Cl ) 



The measurement operator Mq is rewritten as a product 



K 



M Q = -U g M c U g , M e = ]J M Ci (27) 



where the measurement operator 

M Ci =l-2|ci)<Ci| = X E *(l-2|l)„<l| n )A^, (28) 

stabilizes the orthogonal complement of the classical 
state |cj). The components of the binary vector are 
the respective complements of Gj, Cj, = 1 © Ci.j. The 
operator in the parentheses in Eq. (Eq) is the (n — 1)- 
controlled Z gate (the Z-operator is applied to the n-th 
qubit only if all the remaining qubits are in state |1)). 
It can also be represented as n-qubit controlled phase 
gate C!^ = exp(i9P\i\) with 9 = w, where the opera- 
tor P\x) n projects onto the state |1)„. This can be fur- 
ther decomposed as a product of two Hadamard gates 



and an (n — l)-controlled CNOT gate [Fig. g] which for 
71 > 6 can be implemented in terms of one ancilla and 
8(n — 4) three-qubit Toffoli gates |3^] and therefore has 
linear complexity in n. With no ancillas, the complexity 
of the (n - l)-controlled CNOT gate is 0(n 2 ) M. 



ci = 



- Z - 



- H -©- H - 



FIG. 3: Decomposition of n-qubit controlled-Z gate C n in 
terms of (n — l)-controlled CNOT gate (n-qubit Toffoli gate) 
for n = 3. 

The corresponding ancilla-based measurement for M Ci 
can be constructed with the help of two Hadamard gates 
[Fig. ||] by adding an extra control to each gate. In- 
deed, this correlates the state |1) of the ancilla with M Ci 
acting on the n qubits, and the state |0) of the ancilla 
with 1. When constructing the measurement for the 
product of the operators M Ci , it is sufficient to use only 
one ancilla, since for each basis state \wi) only one of 
these operators acts non-trivially. The classical part of 
the overall measurement circuit without the graph state 
encoding Ug is shown in Fig. [|. 

The complexity of measuring Mq [Eq. (p7[)] then be- 
comes K times the complexity of (n + l)-qubit Toffoli 
gate for measuring each M Ci , plus the complexity of the 
encoding circuit Ug and its inverse Ug, which is at most 
n 2 . Overall, for large n, the measurement complexity is 
no larger than n 2 + KO(n), or (1 + K)0{n 2 ) for a circuit 
without additional ancillas. 

We would like to emphasize that so far we have only 
constructed the measurement for error detection. Actual 
error correction for a non-additive code in this scheme 
involves constructing measurements Me = EMqE' for 
all corrupted subspaces corresponding to different degen- 
eracy classes given by different C\g (E) . This relies on the 
orthogonality of the corrupted subspaces, see Appendix 
|A[ For a general i-error correcting code, the number of 
these measurements can reach the same exponential or- 
der B(n, t) as the number of correctable errors in Eq. (0). 
For non-degenerate codes, we cannot do better using this 
method. 

Note that the measurement circuit derived in this sec- 
tion first decodes the quantum information, then per- 
forms the measurement for the classical code, and finally 
re-encodes the quantum state. 



IV. MEASUREMENTS FOR UST CODES 

In this section we construct a quantum circuit for the 
measurement operator Mq of a USt code ((n,K2 k ,d)). 
To this end, we define the logical combinations of non- 
Pauli measurements in agreement with analogous combi- 
nations defined in Ref. p(| for the projection operators, 
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H 
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X° 2 



X CK 



X Ck 



|0> P Mc IV} + \1)Qm c \i>) 
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c 



FIG. 4: Generic measurement of the classical part Mc of the CWS code stabilizer Mq [Eq. (E^)] uses if (n+l)-qubit controlled- 
Z gates. Here X Ci indicates that a single-qubit gate X Ci - j is applied at the j th qubit, j — 1, . . . n, and C Mi is the controlled- M Ci 
gate [see Eq. (p8|)]. This can be further simplified by combining the neighboring Xf^-gates and replacing the controlled- Z gates 
by (n + l)-qubit Toffoli gates as in Fig. til 
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FIG. 5: Measurement Mi A Mo by performing logical AND operation on two ancillas. We use the notations Pi = Pm ( , 
Qi = 1 - Pi, i = 0,1 and assume PiP = PoPi- The circuit returns |,f) = |1) PiP |V>) + |Q) (1 - APo) |^>- Intermediate 
results are: \a) = |000) Pi jV) + 1 100) Qi |^), |6) = |000) PiP |V) + |010) PiQo W + 1 100) QiP \ip) + |110) QiQo |^), |c) = 
j001) PiPo \ip) + 010) PiQo |^> + 1 100} QiPo \ip) + 1 110) QiQo j^}; the last two groups of gates return the first two ancillas to 
the state |00). 
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FIG. 6: Simplified measurement for Mi A Mo; notations as in Fig. Intermediate results (cf. Fig. [l]) are: [o) = ] 00) Pi 1^} + 
1 10} Qi \ip), \b) = |0} (|1> P + j0} Qo)Pi \ip) + 1 10} Qi \ip). The effect of the last block is to select the component with the first 
ancilla in the state |0}; this requires that the projectors Pi and Po commute. The result is equivalent to |1}Pa/ 1 am W) + 

|0) (1 - P Ml AM ) |V}- 



and construct the circuits for logical combinations AND 
[Figs. [| § and XOR [Figs. @, §. We use these circuits to 
construct the measurement for Mq with complexity not 
exceeding 2K(n + l)(n — k). 



A. Algebra of measurements 

Logical AND: Given two commuting measurement op- 
erators Mi and Mo, let Mi A Mo denote the measurement 
operator that stabilizes all states in the subspace 

V(M 1 A M ) = V{Mi) n V(M ). (29) 

The output of the measurement Mi A M is identical to 
the logical AND operation performed on the output of 
measurements M% and Mo. This measurement can be im- 
plemented by the circuit in Fig. |5|. Here the first two an- 
cillas are entangled with the two measurement outcomes; 
the third ancilla is flipped only if both ancillas are in the 
|0) state, which gives the combination |1) PmiAM |V)> 
The projector onto the positive eigenspace of Mi A Mq 



satisfies the identity 

PmiAM = PmiPmq- (30) 

This identity can be used to obtain a simplified circuit 
which only uses two ancillas, see Fig. ^, with the price of 
two additional controlled-Hadamard gates [Fig. 0|. 

The circuits in Figs. || and |6| can be generalized to 
perform the measurement corresponding to the logical 
AND of £ > 1 commuting measurement operators with 
the help of associativity, e.g., M 2 A Mi A M = M 2 A 
(Mi A Mo). The generalization of the simplified circuit 
in Fig. |^ requires only two ancillas for any £ > 1. The 
corresponding complexity is 2£ — l times the complexity 
of a controlled-M gate, plus 2£—l times the complexity of 
a controlled single-qubit gate. When all Mi are n-qubit 
Pauli operators, the overall complexity with two ancillas 
is (2£- l)(n+l). 

Logical XOR: In analogy to the logical "exclusive OR" , 
we define the symmetric difference AAB AC . . . of vector 
spaces A, B, C, ... as the vector space formed by the 
basis vectors that belong to an odd number of the original 
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FIG. 7: Measurement Mi EB Mo by performing logical XOR gate on the two ancillas and subsequent recovery of the first ancilla. 
Notations as in Fig. |. The final result is |/> = |1) (P Qi + PiQo) \i>) + |0> (P Pi + QoQi) \i>)- 
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H 



IV) Mo ~ Afi 



I/) 



FIG. 8: Simplified measurement for Mi EB Mo. Notations as 
in Fig. |. The result (cf. Fig. % \f) = |1) (QxPq + PiQo) IV 1 ) + 
!0) (P1P0+Q1Q0) |V>) is equivalent to |1) P Ml mu \ip) + \0) (t- 



MifflM J 



from that in Fig. 0by a sequence of circuit simplifications 
(not shown). 

The circuit in Fig. ^ is immediately generalized to a 
combination of more than two measurements, Mi EH ... EH 
Mi = (— l) l ~ 1 Mi . . . Mi. The corresponding complexity 
for computing the XOR of t measurements is simply the 
sum of the individual complexities, implying that this 
concatenation has no overhead. 



- H - 



-tt/4 



z 



R 



ir/4 



FIG. 9: Implementation of the controlled- H gate based on 
the identity H = exp(~iYTv/8)Z exp(iY%/8). 



vector spaces. For two vector spaces 

AAB=(An B x ) ffi(Bfl A ± ). 



(31) 



This operation is obviously associative, A A (B A C) — 
(AAB)AC. For two commuting measurement operators 
Mo, Mi, let Mi EH Mq be the measurement operator that 
stabilizes the subspace V{Mi) A V{Mq). Explicitly, 

V(Mi ffl Mo) = V{Mx) A V(M a ) 
= [V(M X ) n V X (M )] © [V(Mo) n V X (M X )] . (32) 

The output of measuring Mi EH Mo is identical to the 
logical XOR operation performed on the outputs of mea- 
surements Mi and Mo. The corresponding measurement 
can be implemented by combining the two ancillas with 
a CNOT gate [Fig. 0. 

To simplify this measurement, we show that Mi ffl 
Af = -MiM . Indeed, Eq. Q implies that for the 
projection operators Pi = Pm^, i = 0, 1, 

P Ml \SM a = - P ) + P {t - Pi)- 

The corresponding measurement operator is factorized 
with the help of the projector identities (|J), 

AfifflAfo = 2[Pi(l-P )+Po(l-Pi)]-l 
= -(2Pi - 1)(2P - 1) = -MiM . 

This implies that V{M X ffl M ) = V X {M\M Q ). In other 
words, the measurement of Mi ffl M can be implemented 
simply as an (inverted) concatenation of two measure- 
ments, see Fig. |8|. The same circuit can also be obtained 



B. Error detection for union stabilizer codes 

A USt code Q = ((n,K2 k ,d)) is a direct sum © of 
K mutually orthogonal subspaces obtained by translat- 
ing the originating stabilizer code Qo = [[^^cfo]]- For 
mutually orthogonal subspaces A _L B, we have A c B 
and B c A , and the direct sum is the same as the 
symmetric difference @), A © B = A A B. 

In turn, the stabilizer code Q is an intersection of the 
subspaces stabilized by the generators Gi of the stabi- 
lizer, see Eq. (0). The translated subspaces tj(Q ) are 
stabilized by the Pauli operators Mjj = 
therefore decompose the USt code Q as 



tjGity We can 



K r-n — k -i K rn—k 

j=l L i=l J j—l L i=l 



(33) 



This gives the decomposition of the measurement opera- 
tor Mq whose positive eigenspace is the code Q as 



M Q = EB 



A M ^ 



(34) 



Recall that the complexity of each of the K logical 
AND operations is [2(n — fc) — l](n + 1). No additional 
overhead is required to form the logical XOR of the re- 
sults. Thus, we obtain the following 

Theorem 1 Error detection for a USt code of length n 
and dimension K2 k , formed by a translation set of size 
K > 0, has complexity at most 2K(n — k){n + 1). 

Note that in the special case of a CWS code (k = 0), 
the prefactor of K is quadratic in n w hereas the corre- 
sponding prefactor obtained in Sec. Ill is linear in n. The 
reason is that in Eq. (p^ ) the graph encoding circuit Qg 
with complexity 0(n 2 ) is used only twice, and the pro- 
jections onto the classical states have linear complexity. 
In Eq. (Ell) we are using K projections onto basis states 



9 



of the quantum code. The advantage of the more com- 
plex measurement constructed in this section is that it 
does not involve having unprotected decoded qubits for 
the entire duration of the measurement. 



V. 



STRUCTURED MEASUREMENT FOR CWS 
CODES 

A. Grouping correctable errors 



Recall from Section II C that for stabilizer codes the 



representatives of the error degeneracy classes form an 
Abelian group whose generators are in one-to-one corre- 
spondence with the generators of the stabilizer. Measur- 
ing the n — k generators of the stabilizer of a stabilizer 
code [[n, k,d\] uniquely identifies the degeneracy class of 
the error. 

In this section we establish a similar structure for CWS 
codes. First, for any subset T> c £ of correctable errors 
of a quantum code Q, we define the set £p of unrelated 
errors which do not fall in the same degeneracy class with 
any error from T>. The formal definition 



% = {E : E e £ A C E ,E X = 0, V£i G V} 



(35) 



is based on the general error correction condition (^) and 
the orthogonality of corrupted spaces, see Appendi x KL 
When errors in £ are non-degenerate, the definition ( pq ) 
is equivalent to the set difference, £\T>. In the general 
case, since we do not distinguish between mutually de- 
generate errors, £p can be thought of as the difference 
between the sets of degeneracy classes in £ and in D. 
Definition (]3q) implies that the subspaces 



V(Q) = E(Q) 



Eev 

and £p(<2), defined analogously, are mutually orthog- 
onal. Moreover, if the elements of the set T> form a 
group 3 = T>, the subspace 3(Q) is also orthogo- 
nal to £-^[3(Q)] [see Eq. (|o|) below]. In other words, 
Qs> = 3{Q) can be viewed as a quantum code which 
detects errors from S-^. 

This observation, together with the error-detection 
measurement for USt codes constructed in the previous 
section, forms the basis of our error grouping technique. 
We prove the following 

Theorem 2 For a CWS code Q — (Q,C) in standard 
form and a group 3 formed by graph images of some 
correctable errors in £, the code Qs> = 3(Q) is a USt 
code which detects all errors in £-^. 

Proof. First, we show that the subspace 3(Q) is a USt 
code. The corresponding set of basis vectors is 

K 

# ({!«*) ,...>*)}) = |J UKW}- (36) 

e Q G® i=l 



These vectors are mutually orthogonal, 

(wi\ e\e a \uij) = 0, Vi,j<K, e a ,ep e 3, e a ^ ep, 

(37) 

since every element e a = Z^ ^^ of the group 3 is a 
representative of a separate error degeneracy class. Fur- 
ther, the group 3 is Abelian, and its elements commute 
with the codeword generators W\ — Zf, Cj € C. There- 
fore, using Eq. (EG), we can rearrange the set (p6[) as 



K 



>{{\w 1 ),...,\w K )}) = {jW l ({e a \s)} 



(38) 



The set in the parentheses on the right hand side is a basis 
of the additive CWS code Qo® formed by the group S> 
acting on the graph state |s). Then, we can write the 
subspace S(Q) explicitly as a USt code [cf. Eq. (p5h] 



K 



#(0 = ©wi(eo»), 



(39) 



where the translations are given by the set of codeword 
operators W = {Wjljl] of the original code Q. Orthog- 
onality condition ( |l4|) is ensured by Eq. (37). 

Second, we check the error-detection condition (Q) for 
the code (|39|). Explicitly, for an error E € and for 
the orthogonal basis states e a Wi \s), 



(s\ Wje^EepWi \s) = ± (wi\Ee\fi P \ Wj ) = 



(40) 



for all a,P,i,j, according to Eqs. (||), ( |35|) and the group 
property of $>. □ 
Now, to correct errors in groups, we just have to find a 
suitable decomposition of the graph images of the original 
error set into a collection of groups, C\g{£) = (J • 31 j, and 
perform individual error-detection measurements for the 
auxiliary codes until the group containing the error 
is identified. 

To find an error within a group 3 = {g\, ... , g m ) with 
m generators, we can try all m subgroups of 3 with one 
generator missing. More specifically, for a generator gi we 
consider the subgroup 3^ = {g\, ... , gr;_i,p/ + i, . . . , g m ) 
and perform error detection for the code Q ( § = 3 {l \Q). 
After completing m measurements, we obtain a repre- 
sentative of the actual error class. This is the product of 
all generators gi for which the corresponding code 
detected an error. 



B. Complexity of a combined measurement 

To actually carry out the discussed program, we need 
to construct the n — m generators Gi of the stabilizer of 
the code Qo@- The generators have to commute with the 
m generators g a in the group 3. 

This can be done with the Gram-Schmidt (GS) orthog- 
onalization p3| of the graph-state generators Si [Eq. Jl6|)] 
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with respect to the generators g a . As a result, we obtain 
the orthogonalized set of independent generators S[ such 



that g a S[ 



-1) dia S' i g a . We can take the last n 



of the obtained generators as the generators of the stabi- 
lizer, G % — S' i+m , i = 1, . . . , n — m. 

The orthogonahzation procedure is guaranteed to pro- 
duce exactly m generators S' a anti-commuting with the 
corresponding errors g a , a = 1, . . . , m. Indeed, the GS 
orthogonahzation procedure can be viewed as a sequence 
of row operations applied to the original n x m binary 
matrix B with the elements bi a which define the original 
commutation relation, 



Sig a = (-l) bic "gaSi 



(41) 



The generator g a ant i- commutes with at least one op- 
erator in .y if and only if the a-th column of B is not 
an all-zero column. Then all m generators are indepen- 
dent (no generator can be expressed as a product of some 
others) if and only if B has full column rank. 

By this explicit construction, the generators Gi of the 
stabilizer of the auxiliary code <2@ [Eq. (|39|)1 are Pauli op- 
erators in the original graph-state basis. The complexity 
of each error-detection measurement M(Q@) is therefore 
given by Theorem |l|. 



C. Additive CWS codes 

The procedure described above appears to be ex- 
tremely tedious, much more complicated than the syn- 
drome measurement for a stabilizer code. However, it 
turns out that for stabilizer codes this is no more diffi- 
cult than the regular syndrome-based error correction. 

Indeed, for a stabilizer code Q = [[n, k,d\], the degen- 
eracy classes for all correctable errors form a group of 
all translations of the code, Si = 3? = (gi, . . . ,g n -k), 
with n — k generators. To locate the error, we just 
have to go over all n - k USt codes @®{Q) gener- 
ated by the subgroups of 9 with the generator gi miss- 
ing. Since the originating code Q is a stabilizer code, 
the USt codes 3>^(Q) are actually stabilizer codes, en- 
coding fc = fc + (n — k — 1) = n — 1 qubits each. 
Hence there is only one non-trivial error, and up to er- 
ror degeneracies, £@(o = (gi)- The corresponding sta- 
bilizers have only one generator each. The neces- 
sary measurements are just independent measurements 
of n — k Pauli operators, the same as needed to measure 
the syndrome. Moreover, if the error representatives gi, 
I = 1, . . . , n — k, are chosen to satisfy the orthogonality 
condition Gigi — {— l) 5 *' g\Gi as in Example 1, the op- 
erators to be measured are the original generators Gi of 
the stabilizer, and the corresponding measurement is just 
the syndrome measurement. 

Example 5. Consider the additive code ((5,2,3)) 
equivalent to the stabilizer code [[5,1,3]], see Examples 
1 and 4. The graph-induced maps of single-qubit er- 
rors form a group of translations of the code, & = 



(Zi, Z 2 , Z 3 , Z4). This group contains all error degener- 
acy classes, S 1 = 2T . With the addition of the logical 
operator X = ZZZZZ, these can generate the entire 
5-qubit Hilbert space Wf 5 from the graph state |s); we 
have ,^(Q) = Uf 5 . 

Indeed, if we form a measurement as for a generic 
CWS code, we first obtain the stabilizer of the auxil- 
iary code Qo@ [Eq. (H^)] which in this case has only one 
generator, 5? = (S5), where S5 = ZIIZX, see Exam- 
ple 4. Translating this code with the set (in this case, 
group) W = {I, X} of codeword operators, we get the 
auxiliary USt code W(Qqs>) as the union of the positive 
eigenspaces of the operators [see Eq. (|32 



Mio = S 5 , M hl =XS 5 X [ 



which is the entire Hilbert space, W{Qqs,) — &{Q) = 
Hf 5 , as expected. 

To locate the error within the group S> with m = 4 gen- 
erators, we form a set of smaller codes = 2#( l \Q), 
I = 1,...,4, where the group Ql^ is obtained from <3) 
by removing the Z-tli generator. The corresponding sta- 
bilizers are = (Si,S 5 ), y {2) = (S 2 ,S 5 ), etc. The 
matrices M^'- of conjugated generators have the form, 



e.g. 



M 



(i) _( S x , -Si 



Sr, 



-s 5 



S2, —5*2 
S5, —S$ 



The code Q% is formed as the union of the common 
positive eigenspaces of the operators in the columns of the 
matrix Mf - . Clearly, these codes can be more compactly 

introduced as positive eigenspaces of the operators G; = 
SiS^, I = 1,...,4. Such a simplification only happens 
when the original code Q is additive. While the operators 
Gi are different from the stabilizer generators in Eq. (11), 
they generate the same stabilizer 5? — (Gi , . . . , G4) of 
the original code Q. It is also easy to check that the same 
procedure gives the original generators Gi [Eq. (pr|)] if we 
start with the error representatives (|13|). O 



D. Generic CWS codes 

Now consider the case of a generic CWS code Q = 
((n,K,d)). Without analyzing the graph structure, it 
is impossible to tell whether there is any set of clas- 
sical images of correctable errors that forms a large 
group. However, since we know its minimum distance, 
we know that the code can correct errors located on 
t = [(d — 1)/2J qubits. All errors located on a given 
set of qubits form a group. Therefore, by taking an in- 
dex set Ac of s < i different qubit positions, 
we can ensure that the corresponding correctable errors 
form a group with 2s independent generators. 
The corresponding graph images C\g(Ej) £ ^4 obey the 
same multiplication table, but they are not necessarily in- 
dependent. As a result, the Abelian group 9a generally 
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has m < 2s generators. Since all group elements corre- 
spond to correctable errors, the conditions of Theorem || 
are satisfied. 

Overall, to locate an error of weight t or less, we need 
to iterate over each (but the last one) of the (™) index sets 
of size t and perform the error-detecting measurements 
in the corresponding USt codes £a{Q) until the index 
set with the error is found. This requires up to (") — 1 
measurements to locate the index set, and the error can 
be identified after additional m < 2t measurements. This 
can be summarized as the following 

Theorem 3 A CWS code of distance d can correct errors 
of weight up to t = \_(d — 1 )/ 2J by performing at most 



N(n,t) 



2t-l 



(42) 



measurements. 



For any length n > 3, this scheme reduces the total num- 
ber (|l|) of error patterns by a factor 



B(n,t) 
N(n,t) 



> 



3n+l 
n+1 

3*, 



, if t = 1, 



if t > 1. 



(43) 



Example 6. Consider the ((5,6,2)) code previously 
discussed in Example 3. While the distance d = 2 
is too small to correct arbitrary errors, we can correct 
an error located at a given qubit. Assume that an er- 
ror may have happened on the second qubit. Then we 
only need to check the index set A = {2}. The errors 
{1, X2, Y2, Z2} located in A form a group with genera- 
tors 1X2,^2}; the corresponding group of classical er- 
ror patterns induced by the ring graph in Fig. ^|(b) is 
$>A = (Z1Z3, Z2). The three generators Gi of the sta- 
bilizer of the originating USt code Qo@ A can be chosen 
as, e.g., Gi = S1S3 = XIXZZ, G 2 = S 4 = IIZXZ, 
G3 = S5 = ZIIZX. Using the classical codewords (|25| ) 
for the translation operators tj — Z Cj , we obtain the 
conjugated generators Mjj = tjGiij 

(Gi, —Gi, Gi, Gi, Gi, — Gi \ 
G 2 , G 2 , -G 2 , -d, G 2 , -d . (44) 
G3, — G3, G3, — G3, — G3, G3 J 

According to Eq. (|33|), the auxiliary code Q@ A is a direct 
sum of the common positive eigenspaces of the operators 
in the six columns of the matrix (fl4|). 

To locate the actual error in this 24-dimensional space, 
we consider the two subgroups = (ZIZII) and 

^( 2 ) = (IZIII) oi&A- The stabilizers of the correspond- 
ing auxiliary codes Q@ and Q 



i( 2 ) 



can be obtained by 

= S 3 = IZXZI, 
respectively; this adds one of the rows 



adding G { 4 1] = S 2 = ZXZII and G^ 



M 



(i) _ 

4,j ^ 



{S2,— S2, 5*2,— ^2, S2,— S2) 



MP- — (S3, — S3,— S3, S3,— S3, S3) 



(45) 
(46) 



to the matrix (Q). The original code Q is the inter- 
section of the codes Q@ an d Q^i t ne corrupted space 
Xi(Q) is located in Q$ , but not in Q^ 1 , while, e.g., 
the corrupted space Y 2 (Q) is located in Q@ A , but not in 

□ 



(i) 



or Qg. 



E. General USt codes 



A similar procedure can be carried over for a general 
USt code ((n, K 2 k , d)), with the only difference that the 
definitions of the groups & and the auxiliary codes Qo® 
[Eq. ({3^) 1 should also include the k ge nerators of the orig- 
inating stabilizer code Qq [Sec. [ID]. Overall, the com- 
plexity of error recovery for a generic USt code can be 
summarized by the following 

Theorem 4 Consider any t-error correcting USt code of 
length n and dimension K2 k , with the translation set of 
size K. Then this code can correct errors using (™) +2t — 
1 or fewer measurements, each of which has complexity 
2K(n + l)(n - k - 1) or less. 



F. Error correction beyond t 

For additive quantum codes, the syndrome measure- 
ment locates all error equivalence classes, not only those 
with "coset leaders" of weight s < d/2. The same could 
be achieved with a series of clustered measurements, by 
first going over all clusters of weight s = t, then s = t+1, 
etc. This ensures that the first located error has the 
smallest weight. In contrast, such a procedure will likely 
fail for a non-additive code where the corrupted spaces 
E\(Q) and Ei(Q) can partially overlap if either E\ or 
Ei is non-correctable. For instance, the measurement in 
Example 6 may destroy the coherent superposition if the 
actual error (e.g., Zi, i ^ 2) was not on the second qubit. 

Therefore, if no error was detected after (") — 1 mea- 
surements, we can continue searching for the higher- 
weight errors only after testing the remaining size-t index 
set. With a non-additive CWS code, generally we have 
to do a separate measurement for each additional cor- 
rectable error of weight s > t. 



VI. CONCLUSIONS 

For generic CWS and USt codes, we constructed a 
structured recovery algorithm which uses a single non- 
Pauli measurement to check for groups of errors located 
on clusters of t qubits. Unfortunately, for a generic CWS 
code with large K and large distance, both the number of 
measurements and the corresponding complexity are ex- 
ponentially large, in spite of the exponential acceleration 
already achieved by the combined measurement. 

To be deployed, error-correction must be comple- 
mented with some fault-tolerant scheme for elementary 
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gates. It is an important open question whether a fault- 
tolerant version of our measurement circuits can be con- 
structed for non-additive CWS codes. It is clear, how- 
ever, that such a procedure would not help for any CWS 
code that needs an exponential number of gates for recov- 
ery. Therefore, the most important question is whether 
this design can be simplified further. 

We first note that the group-based recovery [see The- 
orem D is likely as efficient as it can possibl y be , illus- 
trated by the example of additive codes in Sec. V_C where 
this procedure is shown to be equivalent to syndrome- 
based recovery. Also, while it is possible that for fixed 
K the complexity estimate of Theorem |l| can be reduced 
in terms of n (e.g., by reusing ancillas with measured 
stabilizer values), we think that for a generic code the 
complexity is linear in K. 

However, specific families of CWS codes might be rep- 
resented as unions of just a few stabilizer codes which 
might be mutually equivalent as in Eq. (|l5|), or non- 
equivalent 1 28 . The corresponding measurement com- 



plexity for error detection would then be dramatically re- 
duced. Examples are given by the quantum codes derived 
from the classical non-linear Goethals and Preparata 
codes H||. 

Another possibility is that for particular codes, larger 
sets of correctable errors may form groups. Indeed, we 
saw that for an additive code ((n, 2 k , dj), all error de- 
generacy classes form a large group of size 2™~ fc which 
may include some errors of weight well beyond t. Such 
a group also exists for a CWS code which is a subcode 
of an additive code. There could be interesting families 
of non-additive CWS codes which admit groups of cor- 
rectable errors of size beyond 2 2t . For such a code, the 
number of measurements required for recovery could be 
additionally reduced. 
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implied in the discussion of degenerate errors for CWS 
codes. However, to our knowledge, it was never explic- 
itly discussed for CWS or USt codes. Since our recovery 
algorithms for CWS and USt codes rely heavily on this 
orthogonality, we give here an explicit proof. 

First, consider a stabilizer code Qq. For any Pauli 
operator E = E[E%, there are three possibilities: (i) E 
is proportional to a member of the stabilizer group, E = 
75, where S <E S"o and 7 = i m , m = 0, . . . , 3, (ii) E is 
in the code normalizer J/q but is linearly independent of 
any member of the stabilizer group, and (iii) E is outside 
of the normalizer, E JVq. 

Case (i) implies that the space E(Qq) is identical to the 
code Qo] the errors E\ and Ei are mutually degenerate. 
Indeed, for any basis vector the action of the error 
E \i) = 75 \i) = 7 \i) just introduces a common phase 7; 
any vector € Qq is mapped to 7 \ip) and hence no 
recovery is needed. 

In case (ii) the operator E also maps Qo to itself, but 
no longer identically. Therefore, at least one of the two 
errors E\, E2 is not correcta ble. Indeed, in this case we 
can decompose E (see Sec. II C) as the product of an 



element S £ ,S^o in the stabilizer and logical operators, 
a b 

i.e., E = i m SX Z , where m — 0, ... ,3 determines the 
overall phase. While S € J?q acts trivially on the code, 
the logical operator specified by the binary-vectors a, b 
is non-trivial, wgt(a) + wgt(b) ^ 0. Using the explicit 
basis (0), it is easy to check that the error-correction 
condition (||) is not satisfied for the operators E±, E 2 . 

Finally, in case (iii) the spaces E(Q ) and Q are mu- 
tually orthogonal. Indeed, since E is outside of the code 
normalizer .y^, there is an element of the stabilizer group 
S 6 S^o that does not commute with E. Therefore, for 
any two states in the code, \<p) , \ip) € Qo, we can write 

(<p\ E \4>) = (<p\ ES \i/>) = - (<p\ SE \i/>) =-{tp\E |V) , 

(Al) 

which gives (tp\ E \ip) = 0, and the spaces Qo and E(Qo) 
[also, Ei(Qq) and ^(Qo)] are mutually orthogonal. 

Now, consider the same three cases for a USt code (pf ) 
derived from Q . In case (i) the code is mapped to itself, 
E(Q) = Q. The operator E acts trivially on the code 
(and the errors E\, E % are mutually degenerate) if E 
either commutes (A2) or anti-commutes ([A3]) with the 



entire set of translations generating the code: 



Appendix A: Orthogonality of corrupted spaces 



(Et 3 
or (Etj 



tjE, 



3 = 1, 



-tjE, J = 1, 



,K) 

:,K). 



(A2) 
(A3) 



As discussed in Sec. II B, for a general non-additive 



quantum code Q and two linearly independent cor- 
rectable errors, the corrupted spaces Ei(Q) and E%(Q) 
may be neither identical nor orthogonal [ |17| . However, 
for CWS and USt codes it is almost self-evident that 
when Ei(Q) and E 2 {Q) do not coinside, they are mutu- 
ally orthogonal. This orthogonality is inherited from the 
originating stabilizer code Qq. In particular, in some 
previous publications (e.g., Ref. p0|) orthogonality is 



If neither of these conditions is satisfied, the error- 
correction condition (|^) is violated. This is easily checked 
using the basis = tj \i). 

Similarly, in case (ii), the code is mapped to itself, 
E(Q) = Q, but the error-correction condition (||) cannot 
be satisfied. 

Finally, in case (iii), the space E(Q) is either orthogo- 
nal to Q, or the error correction condition is not satisfied. 
The latter is true if E is proportional to an element in 
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one of the cosets t^tj'^o, where j ^ j', 1 < < K. 
Then the inner product (j, i\ E \j', i) ^ 0, i = l,...,fc, 
which contradicts the error-correction condition (|^). In 
the other case, namely, when E is linearly independent of 
any operator of the form tjtj>S, j,j' = 1, . . . , K, S E S^q, 
E must be a member of a different coset t a 5f§ of the 
stabilizer S^q of the code Qo in 3P n . This implies orthog- 
onality: 

E \f,i') = (i\ t]Et r |0 = i™ (i\ t]t r t a |0 = 0, 



where m — 0, ... ,3 accounts for a possible phase factor. 



Overall, as long as the error correction condition (||) 
is valid for a USt code Q and the Pauli operators E\, 
Ei, the spaces E\{Q) and i?2(Q) either coincide, or are 
orthogonal. Since CWS codes can be regarded as USt 
codes ori ginat ing from a one-dimensional stabilizer code 
Qo, [Sec. HF|, the same is also true for any CWS code. 
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